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ABSTRACT 

We have performed deep X-ray observations of the remnant of Kepler's supernova 
(SN 1604) as a Key Project of the Suzaku Observatory. Our main goal is to detect 
secondary Fe-peak elements in the SN ejecta to gain insights into the Type la supernova 
explosion mechanism and the nature of the progenitor. Here we report our initial results. 
We made a conclusive detection of X-ray emission lines from highly ionized Mn, Cr, 
and Ni as well as Fe. The observed Mn-to-Cr line flux ratio is ~0.60, ~30% larger 
than that measured in Tycho's remnant. We estimate a Mn-to-Cr mass ratio of ~0.77, 
which is strongly suggestive of a large neutron excess in the progenitor star before the 
onset of the thermonuclear runaway. The observed Ni-to-Fe line flux ratio (~0.03) 
corresponds to a mass ratio of ~0.06, which is generally consistent with the products of 
explosive Si-burning regime in Type la explosion models, and rules out contamination 
from the products of neutron-rich nuclear statistical equilibrium in the shocked ejecta. 
Together with the previously suggested luminous nature of the explosion, these mass 
ratios provide strong evidence for a super-solar metallicity in the SN progenitor (~3Z ). 
Kepler's supernova was likely the thermonuclear explosion of a white dwarf formed in 
the recent past that must have exploded through a relatively prompt channel. 
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Subject headings: ISM: supernova remnants — ISM: individual objects (Kepler's super- 
nova remnant) — X-rays: ISM 



1. INTRODUCTION 



The progenitors of Type la supernovae (SNe) are believed to be C/O white dwarfs (WDs) 
in close binary systems. The WD becomes unstable when it approaches the Chandrasekhar limit 
{Mch = 1-44 Mq) by either accreting a sufficient amount of material from a non-degenerate com- 
panion {single- degenerate channel) or merging with another WD {double- degenerate channel). This 
mass increase triggers a thermonuclear runaway and an explosion ensues by the fusion of carbon 
and oxygen. While this is a generally accepted picture, the nature of the progenitor and its binary 
companion and the details of the explosion process remain elusive. The progenitor's nature (e.g., 
age, mass, and metallicity) , its environment (e.g., ambient structure and the companion's nature), 
and the details of the explosion mechanism (whether the burning front is a deflagration, a delayed 
detonations, etc.) may significantly vary among SNe la. Understanding the diversity of SN la is 
essential to reduce systematic uncertainties in their cosmological applications. It has been estab- 
lished that SN la progenitors have a wide range of delay-times, from a few hundred Myr to several 
Gyr, and that "prompt" SNe l a produce more 56 Ni and are brighter than the "delayed" population 
(e.g., iMaoz Mannucci 20121 ). If SN la exploded with delay times of only a few hundred Myr, 
their progenitor's ma sses must have been significantly larger (~3.5-8 Mq) than the Sun in their 
main sequence stage ([Aubourg et al. 2008 ). 



The progenitor's metallicity {Z) is a key parameter related to the age population of SNe la. 
During the main sequence evolution of an intermediate mass star, the CNO-cycle converts C, N, 
and O into 14 N, which subsequently produces 22 Ne through the intermediate (3 + decay of 18 F in 
the star's hydrostatic He-fusion phase. This results in neutron excess (r/ = l-2Y e where Y e is the 
proton to nucleon number ratio) in the progenitor. Except for the deepest core of M < O.2M 
where 77 is dominated by ele ctron-capture in neutro n-rich nuclear statist ical equilibrium (n-rich 
NSE) ([Brachwitz et al. 20 00). 7/ is proportional to Z (jTimmes et al. 20031 ) . During the explosion, 
the excess amount of neutrons is stored in trace elements with unequal number of protons and 
neutrons. 55 Mn is the most abundant among such elements. 52 Cr is produced by the same explosive 
Si-burning as Mn, but its sy nthesis is unrelated to p. Thus, the Mn-to-Cr mass ratio is an excellent 
tracer of the progenitor's Z (jBadenes et al. 20081 : iBravo 20131 ) . The progenitor's Z, together with 
the details of the explosion physics, can a lso influence the synthesis of stable Ni (mostly 58 Ni, e.g., 
Iwamoto et al. 19991 : iTimmes et al. 20031 ). Observational constraints on the content of these Fe- 
group ejecta elements in SNe la would thus provide insights on the detailed nature of the progenitor 
system and the explosio n mechanism. Howey er, the relatively long decay-times of radioactive Fe- 
peak elements (~yr, e.g.. iKuchner et al. 1994 ) make it difficult to detect these elements directly in 
SNe la. Instead young supernova remnants (SNRs) provide a useful opportunity. 



Kepler's SNR ("Kepler" hereafter) is the remnant of SN 1604. The historical SN light curve sug- 
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gested a Type la explosion (|Baade 19431 ). although this was not a unique classification (jDoggett Branch 19851 ) 
Kepler shows strong atomic X-ray emission lines from h i ghly-ionized ions of Si, S, Ca , and Fe (e.g., 
Decourchelle Ballet 1994! ; iKinugasa Tsunemi 1999 : ICassam-Chenai et al. 2004! ). These over- 
all X-ray spectral characteristics are sim ilar to those of Tycho's SNR, an unambiguous Type la 
(jKrause et al. 20081 : iBadenes et al. 20 06). Reynolds et al. (2007) found that the X-ray ejecta are 
dominated by Si, S, and Fe, without evidence for O, Ne, and Mg. Bas ed on this e jecta abundance 
structure and the absence of a conspicuous stellar remnant (see also iBlair 20051 ) , they ruled out 
a core-collapse origin for Keple r. Kepler appears to be interacting with a modified circumstellar 
medium (CSM. iBandiera 19871 ). Overabundant nitrogen was detected in Kepler, indicatin g the 
presence of dense CSM produced by massive stellar winds (jDennefeld 19821 : IBlair et al. 199ll ). The 
spatial distribution o f interacting CSM suggested that the stron g winds might have originated from 
an AGB companion (jChiotellis et al. 20121 : iBurkev et al. 20131 ). These results generally support a 
single-degenerate origin for Kepler. 

Based on observations at seve ral wavebands, a range of distances (from d ~ 3 kpc to d > 6 kpc) 
to Kepler h as been suggested (e.g.. lReynoso fe Goss 19991 : ISankrit et al. 20051 : 1 Aharonian et al. 20081 : 
Vink 20081 ). Recently, Patnaude et al. (2012) showed that hydrodynamic and X-ray spectral simu- 
lations are consistent with the Chandra spectrum of Kepler for d ^> 7 kpc if the shock is expanding 
into massive stellar winds, while the distance may be closer [d ~ 4 kpc) for the case of a uniform 
medium. The Type la explosio n must have been brig ht and 56 Ni-rich to explain the strong Fe line 
emission observed by Chandra (jPatnaude et al. 20121 ). 



With its Type la origin, young age (409 yr), and the ejecta-dominated spectrum of X-ray 
emission, Kepler provides an excellent opportunity to study the details of Type la explosion mech- 
anisms and the progenitor's nature. In this Letter, we report th e initial results fro m our Suzaku Key 
Project of Kepler. As we presented in our preliminary report (jPark et al. 20121 ). we detect X-ray 
emission from highly-ionized Mn, Cr, Ni, and Fe ions. Line flux ratios of these Fe-group elements 
strongly suggest a super-solar Z for the progenitor. We describe our observations in Section [2j Our 
data analysis and results are presented in Section [3l A discussion and conclusions are presented in 
Sections HI 



2. OBSERVATIONS 

We performed our Suzaku observations of Kepler as part of a Key Project in 2009 August-2011 
March during Cycles 4-5 (Table [T{0. The X-ray Imaging Spectrometer (XIS) was chosen to be the 
primary detector. The FI-CCDs (XIS0 and 3) have a better energy resolution (FWHM ~ 180 eV at 
E ~ 6 keV) than the BI-CCD (XIS1, FWHM ~ 250 eV at E ~ 6 keV). The non-X-ray background 
flux is significantly lower (by a factor of ~2-3 at E = 5-8 keV) in the FI-CCDs than the BI-CCD. 



We also used our earlier Suzaku data (ObsID 502078010), performed on 2008 February. 
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When the two available FI-CCDs (XISO and 3) are combined, the collecting area is larger by a 
factor of 2 than that of XIS1. Thanks to the good energy resolution and low background, the FI- 
CCDs efficiently detect and measure faint line emission fluxes from Fe-group ejecta elements at E 
~ 5-8 keV, which is our primary goal. Thus, we used the XISO and 3 data in this work. Some time 
intervals in observations taken in March 2011 were excluded from the analysis because of temporary 
on-board software issues with the XISO (for ~25 ks in ObsID 505092050). Otherwise, we followed 
the standard data reduction process including charge injection, and combined all 3x3 and 5x5 
mode spectra for each XISO and XIS3. After the data reduction, the total effective exposures are 
~665 and ~691 ks for XISO and XIS3, respectively. 

Since Kepler is a bright extended source (~2' in radius), scattered photons from the SNR 
dominate nearly the entire CCD due to the large point spread function of the telescope, which 
makes it difficult to adequately estimate the background spectrum from surrounding regions. The 
characterization of the background spectrum is critical to accurately measure fluxes from faint 
line features. To mitigate the bright scattered light from the SNR, we have performed separate 
background observations. We chose four nearby source-free regions (~1.5° from Kepler), one each 
of which is to north, south, east, and west of Kepler (Table [1]). We processed these background 
observations in the same way as we did for the source observations. After the data reduction, the 
total effective exposure for the background observations is ^273 ks each for XISO and XIS3. We note 
that recently Yang et al. (2013) analyzed our Suzaku data of Kepler using background estimates 
from source- free regions on the same XIS chip instead of using separate background pointings. Their 
measurements of lower fluxes (by ~40%) for Mn Ka and Cr Ka lines than ours (see Section [3]) are 
likely caused by the scattered X-ray contamination in their background estimates. They suggested 
a carbon-deflagration for Kepler's explosion based on a low value of the Cr-to-Fe Ka line equivalent 
width (EW) ratio (~0.008). Our larger Cr line flux indicates the Cr-to-Fe Ka line EW ratio of 
~0.016 which would rather favor delayed detonations for Kepler as are self-consistently suggested 
by our Ni-to-Fe mass ratio measurement (Section 2]). 



3. ANALYSIS & RESULTS 

We extracted the source spectrum from the entire SNR (a circular region of 4' in radius). We 
extracted the background spectrum from the same detector region as the source in each background 
pointing, and combined them. In the 5-8 keV band, we obtained ~130000 counts from the source 
for each of XISO and XIS3. The background contribution is ~7%. For the spectral model fits, 
we re-binned the source spectrum to contain a minimum of 50 counts per energy channel. Our 
main goal is to measure line fluxes from Cr, Mn, Fe, and Ni using the integrated X-ray spectrum 
of Kepler. Thus, we fit the background-subtracted spectrum with a phenomenological model in 
the 5-8 keV band. We fit the observed spectrum with an absorbed power law (PL) model with 
5 Gaussian components. The PL component is to characterize the underlying continuum. Each 
Gaussian is to fit the emission feature from Cr Ka (E ~ 5.5 keV), Mn Ka (E ~ 5.9 keV), Fe 
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Ka (E ~ 6.5 keV), Fe K/3 (E ~ 7.1 keV), and Ni Ka (E ~ 7.5 keV) lines, respectively. The 
same Gaussian line width was assumed for all five lines , and was varied to g ive the best fit. The 



foreground column is fixed at Nr = 5.2 x 10 21 cm -2 ([Reynolds et al. 20071 ). The photon index 
(r) and normalization parameters for the PL and Gaussian components are varied freely. We fit 
the XISO and XIS3 spectra simultaneously, tying all parameters between the two spectra. 

The observed spectrum and the best-fit model (T = 2.60, X / u = 1469.8/1513) are shown in 
Figure [TJ In the 5-8 keV band, X-ray emission is dominated by the strong Fe Ka line at E ~ 
6.45 keV. X-ray emission lines from He-like ions of Mn and Cr are clearly detected (at ~9<7 and 
~14cj level for Mn and Cr, respectively) (Table [2] & Figure [JJ. Fe K/3 (E ~ 7.1 keV) and Ni 
Ka (E ~ 7.5 keV) lines are also conclusively detected. To verify that the Ni line was not due to 
the instrumental Ni fluorescent line, which varies with the incidental cosmic-ray flux, we repeated 
our Ni line flux measurement using an off-source background region from the on-target data. We 
found that the Ni line was clearly present with this background subtraction. We also compared the 
instrumental Au Lya line flux (E = 9.671 keV) between the source and background observations. 
These tests indicated only a ~ 10-30% difference in the instrumental line fluxes between the source 
and background pointings. Thus, we conclude that the detected Ni line originates from Kepler. 
The small difference in the instrumental background does not affect our discussion in Section SJ 
Alternatively, we fit the underlying continuum with a thermal bremsstrahlung model (kT = 5.11 
keV, x 2 / u = 1465.2/1513). The estimated fluxes for Fe and Ni lines are identical to those from the 
PL model fit. Line fluxes from Mn and Cr are estimated to be slightly smaller (by ~10%), but the 
flux ratio is consistent (within la uncertainties) with that estimated by the PL model fit. 

The archival Chandra data of Kepler show thin non-thermal filaments along the outermost 
boundary, which are identified as synchrotron radiation from the shock-accelerated relativistic 



electrons ([Reynolds et al. 20071 ). To test the effect from synchrotron emission in the measured line 
fluxes, we repeated our spectral model fit with a PL + thermal bremsstrahlung + 5 Gaussians. 
Based on PL model fits of several non-thermal filaments using the archival Chandra data of Kepler, 
we estimated T ~ 2.4-2.6 for these features. Thus, we performed the model fits with T fixed at 2.4, 
2.5, and 2.6. Using the archival Chandra data of Kepler, we estimated that the thin non-thermal 
shell (typically with a width of ~5") all around the SNR contributes ~10% of photon counts in 
the 5-8 keV band. Thus, we assumed a 10% fractional contribution from the PL component in 
the 5-8 keV band total flux. These model fits are statistically good (\ 2 /v ~ 1.0) with negligible 
changes in the best-fit parameters from those described above. The best-fit thermal bremsstrahlung 
temperature is kT = 5.1 keV, and the estimated line fluxes are consistent with the results derived 
above. We conclude that our measurements of individual line fluxes are robust, and hereafter our 
discussion is based on the results listed in Table [2j 
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DISCUSSION & CONCLUSIONS 



Our deep Suzaku XIS data allow us to accurately measure the Fe Ko line center (Table 2), 
which constrains the ionization timescale to n e t ~ 2 (±0.03) x 10 10 cnr 3 s for temperatures in the 
range kT = 3-8 keV for the Fe-zone ejecta gas. The estimated Mn-to-Cr flux ratio is /Mn/Zo = 
0.60=fc0.16 (2a uncertainties, hereafter). This flux ratio is ~30% larger than that inferred for Tycho 
(jBadenes et al. 20081 : iTamagawa et al. 20091 ). The Mn-to-Cr mass ratio can be calculated using 
their line flux ratio: Mm d / Mq x = 1.057 (/Mn//cr)/(eMn/£Cr)) where eMn/eo is the ratio of specifi c 
emissivities per ion, and 1.057 is the atomic mass ratio between Mn and Cr (jBadenes et al. 20 08). 



We calculate £Mn/ecr = 0.82±0.20 appropriate for the plasma with kT = 3-8 keV and n e t = 2 x 
10 10 cm -3 s using the new atomic da ta (Eriksen et al. in preparation) generated for this project 



with Flexible Atomic Code (|Gu 20081 ). The uncertainty of emu/^Ct depends on the measurements 
on the gas temperature, ionization timescale, and the atomic data. Based on our inferred ranges 
of kT and n e t as well as the uncertainties of the atomic data, we estimate ~20-25% uncertainties 
on SMn/ecr- Thus, we adopted a conservative limit of a 25% uncertainty on £Mn/ecr- Then, we 
estimate M Mn /M C r = 0.77 



+0.53 
-0.31- 



Assuming that the neutron excess in Kepler's progenitor was not significantly affected by 
the innermost NS E region, we may es timate the progenitor's Z using a PL relation betweeen 
Myi n / Mq t an d Z (IBadenes et al. 20081 ): i.e., unless the non-explosive C-burning, the so-called 
C- simmering (|Piro Bildsten 20081 ). was widespread shortly be fore the SN explosion, Z shows a 
simple relationship of Mmu/ Mq t = 5.3Z 65 (jBadenes et al. 20081 ). While C- simmering may modify 
the neutro n excess in sub-lumin ous SNe, where the explosive Si-burning region extends deeper into 
the ejecta (jBadenes et al. 20081 ). there is no evidence for a significantly sub -luminous explosion for 
Kepler in the historic al light curve. Recent distance estimates of d ^ 6 kpc (jAharonian et al. 20081 : 
Patnaude et al. 20121 ) support a normal or luminous SN for Kepler. The C handra X-ray spectrum 
of Kepler is consistent with an explosion model that produced ~1 M Q of Fe (jPatnaude et al. 20121 ) . 
Thus, Kepler is unlikely the remnant of a sub-lu minous SN. Then, we estimate that the progenitor's 
metallicity is Z/Z^ = 3.6 +tn for Z Q = 0.014 (jAsolund et al. 2009h . and 2.7± 3 J for Z Q = 0.019 



(jAnders &: Grevesse 19891 ) . This result provides the first observational evidence that the progenitor 
of a SN la was a star with a super-solar Z (Figure [2^) . 

We detect the Ni Ka line for the first time in Kepler. Since the bulk of Fe is created by 
the radioactive decay of 56 Ni, M^/M-p e ~ M5s Ni /Ms6 Ni = 1.051 (/Ni//Fe)/(^Ni/eFe), where 1.051 
is the atomic mass ratio between 58 Ni and 56 Fe. For the measured fm/ f~Fe = 0.031±0.005 and 
ENi/eFc = 0.51±0.13 for the Ka line in the same plasma condition discussed above, we estimate 
M 5 8 Ni /M 5 6 Ni ~ 0.06±g;g|. This Ni mass ratio is in good agreement (within ~10%) with delayed 
detonation models, w hile being significantly different (by a factor of >2) from deflagration models 
(jlwamoto et al. 19991 ). 



Our measured M^/M-pe indicates that the bulk of the shocked Fe-group ejecta was synthesized 
in the explosive Si-burning regime, especially considering the high Z inferred from Myi n / Mq t 
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(Figure[2b). Contamination by a small amount of NSE material is possible, but a large contribution 
of products from the n-rich NSE regime deep in the core of the exploding star is ruled out, regardless 
of the Z value. This strengthens our confidence on Myi n / Mq t as a tracer of the neutron excess in 
the progenitor (and hence Z), since it is unlikely that material synthesized in the inner n-rich core 
can contaminate the shocked ejecta and increase My[ n /MQ r while keeping M^i/Mp e at its observed 
value. 

The observed Ka line flux ratio fm/ fFe is higher (by ~25%) in the northern half of Kepler 
than in the southern half, while it is identical (within ~2%) between the eastern and western halves 
of the SNR. The Ni Ka-to-Fe K/3 line flux ratio is also ~70% higher in the northern half than in 
the south, while it is fully consistent between the east and west (within ~2%). The presence of the 
bright bow shock-like X-ray emission feature in the northern shell of Kepler suggests a more efficient 
development of the reverse-shock there than in the southern shell. The reverse-shock might have 
progressed further into the SNR in the northern shell (than in the south) to shock more Ni-rich 
ejecta close to the SNR center. Thus, although the spatial variation of the measured Ni-to-Fe line 
flux ratios is statistically marginal (~1.5-2.5ct), it would not be surprising if such a variation has 
a real physical origin. On the other hand, we do not find evidence for a similar spatial variation 
in the Mn-to-Cr line flux ratio. The Mn-to-Cr line flux ratio is consistent within lcr uncertainties 
between the north/south and the east/west halves of the SNR. This suggests that the correlation 
between the observed Mn-rich and Ni-rich ejecta is probably insignificant even if the Ni-rich ejecta 
were produced in the NSE. Thus, the bulk of the detected neutron excess appears to be associated 
with material synthesized in the explosive Si-burning regime of a relatively bright SN la, and it 
should be a tracer of the progenitor's Z. 

The implied high Z for the progenitor is consistent with the location of Kepler close to the 
Galactic center, where the interstellar chemical composition is expected to fu rther enriched than in 



the solar neighborhood. Since higher- Z stars likely result in a lower-mass core (IWeiss Ferguson 20091 ). 
Kepler's progenitor WD probably went through an efficient mass growth before exploding as a rela- 
tively young prompt population SN la. A quantitative assessment of the suggested spatial variation 
of the Ni-to-Fe mass ratio and the effect from the core electron-capture on the observed neutron 
excess may require detailed hydrodynamic calculations beyond the scope of this Letter. 
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Table 1. Suzaku Observation Log of Kepler's Supernova Remnant 
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262.6730, 


-21 


,4605 


146.2 


505092050 


2011-3-8 


Kepler SNR 


262.6732, 


-21. 


.4607 


84.6 


505092060 


2011-3-14 


Kepler SNR 


262.6734, 


-21. 


,4607 


46.5 


505092070 


2011-3-29 


Kepler SNR 


262.6738, 


-21. 


.4631 


133.4 



Table 2. Summary of Detected Atomic Emission Lines 



Line Center Energy 


Line Width a 


Line Flux 


Element 


(keV) 


(eV) 


(10~ 6 photon cm" 2 s" 1 ) 




5.50±0.02 


65.5±1.2 


8.31±1.20 


Cr (Ka line) 


5.99±0.04 


65.5±1.2 


4.99±1.16 


Mn (Ka line) 


6.449±0.0008 


65.5±1.2 


365.91±2.78 


Fe (Ka line) 


7.12±0.02 


65.5±1.2 


9.92±1.24 


Fe (K/3 line) 


7.53±0.02 


65.5±1.2 


11.18±1.88 


Ni (Ka line) 



Note. — 2a uncertainties are shown. 
a The line width is tied among all lines. 
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Fig. 1. — Suzaku XIS spectrum of Kepler's SNR. The best-fit model (PL + 5 Gaussians) are 
overlaid. 
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Fig. 2. — (a) My \ n /Mc,r as a function of Z in SN la explosion models (adapted from 
Badenes et al. 200811. The mod els include ID and 3D simulations for all major explosion mecha- 



nisms (see lBadenes et al. 20081 ). The range of solar Z and the corresponding Mmu/ Mq v are shown 
by the grey area. Our measured Mmh/Mct for Kepler is marked. For comparison, Myin/Mci for 
Tycho's SNR is also shown to the right of the plot, (b) M^/M-p e as a function of progenitor Z for 
different burning re gimes for delayed detonation SN la explosion models (the same models from 



Badenes et al. 20081 ). Our estimate for Kepler is shown. 



